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Abstract. We explain how the kinetic theory of L. Boltzmann is applied to weakly nonlinear wave 
equations. 



Introduction 

Kinetic Theory is an everlasting contribution of Ludwig Boltzmann to physics. While 
he devised the theory to understand the dynamics of dilute gases, the method as such is 
much more general and has been applied with great success in many areas. A more recent 
example is the flow of granular media. Their defining property are inelastic collisions. 
Thus a natural, in fact highly informative, approach is to follow the route of kinetic theory, 
adjusting for the appropriate collision mechanism (TJ. 

In the early days of kinetic theory it was a convoluted process to understand that the 
Boltzmann equation cannot possibly hold for every mechanical initial condition. But to 
set forward the proposition that the Boltzmann equation becomes exact in a particular lim- 
iting procedure, now called the Boltzmann-Grad limit [2], took almost 80 years counted 
from Boltzmann's 1872 paper. In this limit the sometimes painfully vague notions of 
"overwhelming probability" and the like acquire a definite meaning. For example, the 
set of exceptional initial phase points has a measure (w.r.t. the initial measure on phase 
space) which tends to zero in the Boltzmann-Grad limit. Grad's proposition was proved 
as a mathematical theorem by Lanford |3, see 0121 for more expanded discussions. 

Given the example of low density gases, one might wonder whether similar ideas are 
applicable to further microscopic systems. The Landau equation for weakly interacting 
particles and the Vlasov equation for weak, long range forces are examples from plasma 
physics. After the advent of quantum mechanics, an obvious task was to also adjust 
kinetic theory. For electrons, satisfying Fermi-Dirac statistics, the generalization was 
accomplished by Nordheim [6 1, for the quantized wave equation (phonons) by Peierls Q, 
and for dilute quantum gases by Uehling and Uhlenbeck (8). 

The derivation of the Boltzmann equation for interacting quantum particles is well 
covered in recent articles Il9l flOl [TTII . Here I summarize the current status for weakly 
nonlinear wave equations with random initial data. On the surface this problem looks 
simpler than weakly interacting quantum particles, since there is no need for operators, 
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Fock spaces, and the like. On the other hand, at least on the level of time-dependent 
perturbation theory, the wave equation and quantum systems have similar Feynman dia- 
grams. In particular cases the wave diagrams are a proper subset of the set of all quantum 
diagrams. As a consequence the limiting kinetic equations differ only minimally. For 
the kinetic description a dividing line seems to be point particles versus wave equation, 
the latter case being much less understood because of the nonlocal interaction once the 
dynamics is expressed in terms of multi-point Wigner functions. 

Writing down a kinetic equation always involves some modeling aspects. The avail- 
able experience provides sufficient guidance as to what would constitute an acceptable 
kinetic equation. There is then no point in waiting for a firm microscopic derivation. 
Rather more urgently is to work out predictions from the kinetic equation, which by it- 
self is a demanding endeavour. In fact, a kinetic equation to be exact in a particular limit 
might very well be exceptional. A prominent example are lattice Boltzmann algorithms, 
which are used widely to simulate flows in complicated geometries and close to react- 
ing surfaces. For such algorithms the Boltzmann equation remains as an approximation, 
although a very powerful one. 

To give a brief outline, in the following section we explain, what I call, the kinetic 
framework. Physicists tend to be more computationally oriented and the underlying sta- 
tistical properties of the microscopic system are rarely spelled out in full detail. But on 
a conceptual level the kinetic framework is most useful and, assuming it to hold, the ki- 
netic equation is obtained rather easily. We emphasize that the kinetic framework is not 
restricted to a particular model. It is a general structure, which can be used whenever 
a splitting into a "free" part and an appropriately small perturbation is meaningful. In 
Section [2] we write down the phonon Boltzmann equation in case of a quartic potential. 
Some of the basic properties of this equation are discussed in Section [3] Perhaps some- 
what unexpected, the kinetic approach yields interesting dynamical information even for 
an anharmonic chain. In this light we revisit the famous Fermi-Pasta-Ulam problem. 



1 The kinetic framework 

We plan to study the scalar wave equation in three space dimensions with a cubic nonlin- 
earity of the form 

d 2 

— u = (A - ul)u - yflu 3 . (1.1) 
of 

Here u : R x R 3 — » R is the wave field, A the Laplacian, and A > 0, but small. The initial 
data are random and of finite energy. As will be explained, the total energy is taken to 
be of order A~ 3 and the support of u of linear extent /T 1 . To avoid the issue of ultraviolet 
divergencies, Eq. dl.lt is discretized spatially through the standard grid Z 3 . It is important 
to also allow for a general dispersion relation of the linear part. Physically, the discretized 
wave equation describes the lattice dynamics of dielectric crystals, see Ifl2l [131 [141 for 
systematic expositions. In this context the scalar equation would correspond to a one- 
band approximation. 
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Fourier transform will be a convenient tool. Let T 3 = [-A, i] 3 be the first Brillouin 
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zone of the lattice dual to Z . For / : Z — > R its Fourier transform, /, is defined by 

f(k)=Y j e- tMx f x . (1.2) 

.teZ 3 

Here k e T 3 and / extends periodically to a function on R 3 . The inverse Fourier transform 
is given by 

f x = f dke a " kx f{k) . (1.3) 
Jt 3 

For x e Z 3 the displacement away from x is denoted by g A - e R and the corresponding 
momentum by e R. In the harmonic approximation the interaction potential is 



Uhmmiq) = ^ ^ a(x - y)q x q y . (1.4) 
The elastic constants a(x) satisfy 
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a(x) = a(-x) and |a(jc)| < y <?~ riW (1.5) 

for suitable jq, j\ > 0. Mechanical stability requires 

a(k)>0. (1.6) 

In fact, we assume 'a(k) > for k + 0. The case (?(£) A; 2 for small is refered to as 
'acoustical' while c?(0) > is 'optical'. The dispersion relation of the harmonic part is 

co(k) = yfa(kj. (1.7) 

Following dl.lt as a concrete example, we consider a quartic on-site potential. Then 
the full hamiltonian of the lattice dynamics is given by 

H = \ Z & + \ Z a(x - + \ ^ Z ^ ■ (L8) 

xeZ 3 x.yeZ 3 jteZ 3 

The harmonic part of // will be denoted by H^, H - H^ a + y~AV. In a sense to be 
explained below, yfAV is a "small" perturbation of H^ d . 

We concatenate q x and r into a single complex-valued field a{k) as 

fl(fc) - 4=( yf^kJqXk) + i-±—pW) (1.9) 
V2 V V^(fc) 7 

with the inverse 

q(k) = + «(-*)*) . = -r= V^W( - <k) + a(-k)*) . (1.10) 

V2 W V2 
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Then the hamiltonian reads 

H = H h . d +VAV, ff ha = \ dk(o(k)a(k)*a(k) , (1.11) 

Jt 3 



V = — I dk\dkidk-idk\5{k\ + /q + £3 + fct) 
4 

X ]~[ {2cj{kj)Y ll2 {a(kj) + a(-kj)*) . (1.12) 

7=1 

Hence the equations of motion for the a-field are 

— a(k, t) — —ia>(k)a(k, t) — i V/i I dk\dk2dk-i,5{k — k\ — ki — k-i) 
at J T 9 

3 

x(2aj(k)y 1/2 Y\ (2io(kj)y 1/2 (a(kj, t) + a(-kj, t)*) . (1.13) 
7=1 

In particular for A - 0, 

a(k, t) = e~ iajm a(k) . (1.14) 

Let us now briefly recall the case of dilute gases. The central quantity is the Boltzmann 
/-function, which is the number density on one-particle phase space R 3 x R 3 . It changes 
in time by the free motion of particles and through pair collisions, 

1 rft + V-V r ft = C(j t ,f t ). (1.15) 

We do not write out the collision operator explicitly but note that it is bilinear in f, and 
strictly local in the space-variable r, since on the kinetic scale particles collide at the same 
point. At low density the full particle statistics is close to Poisson, jointly in positions 
and velocities, with an intensity given again by the /-function. This double meaning of 
the /-function is a source of conceptual confusion. In the StoBzahlansatz one assumes 
independent incoming velocities, thus the second meaning of /, while the predictions of 
kinetic theory are based on a law of large numbers, thus the first meaning of /. The 
Poisson statistics appears naturally through the good statistical mixing properties of the 
collisionless dynamics. In fact, at infinite volume, these are the only translation invariant 
measures, which are invariant under the collisionless dynamics and which have a finite 
number, energy, and entropy per unit volume lfT31 . 

For wave equations the Poisson measure will be substituted by a Gaussian measure, 
constrained to be locally invariant under the dynamics generated by Hh d . In the spatially 
homogeneous case, imposing time-stationarity and using the explicit solution (11.14b for 
A — 0, the {a(k), k e T 3 } are then jointly Gaussian with mean zero and covariance 



(a(k)*a(k')) = W(k)5(k - k') , 



(1.16) 
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(a(J0*a(kf)*) = , (a(k)a(k')) = , (1.17) 

defining the power spectrum W(k). Clearly, W(k) > 0. 

For the spatially inhomogeneous case the construction is slightly more elaborate. 
Since the anharmonic potential is of order V^, the mean free path for phonons is of 
order A~ l . This is the scale on which spatial inhomogeneities have to be imposed. Local 
stationarity still implies mean zero and ( 11.17b . For the covariance ( I1.16l l we first prescribe 
the limiting local power spectrum W : R 3 x T 3 — > R with W > and a rapid decay in 
position space. In analogy to the semiclassical limit of the Schrodinger equation we refer 
to W as Wigner function. A locally stationary Gaussian measure, denoted by (•) ' , is 
in fact a scale of probability measures depending on A. For each A the measure (-) G,A is 
Gaussian with 

{a(k)f' A = , {a(k)a(k')f' A = . (1.18) 
Following Wigner the local power spectrum of the a-field is defined through 

P\x,k) = 2- 3 \ drie anxn a{k-ril2ya(k + ril2). (1.19) 

J(2T) 3 

We rescale the lattice to have lattice spacing A through the substitution x = A _1 y, y e 
(/1Z/2) 3 , and obtain the rescaled local power spectrum 

P A (y, k) = (A/2) 3 \ drje i2ny il a(k - Aj]/2)*a(k + Arj/2) . (1.20) 

J(2T/^) 3 

P A is a random field. By definition, its average is the one-particle Wigner function, 

W A (y,k) = (P A (y,k)f- A (1.21) 
and its variance the two-particle Wigner function, 

W A (y u k u y 2 ,k2) = (P A (yuk l )P A (y2,k 2 )) c - A . (1.22) 

With L-Jyi denoting modulo A, the scale of Gaussian measures (•) G-/1 is assumed to satisfy 
the pointwise limit 

lim W A (lr} A , k) = W(r, k) . (1.23) 

A — >0 

Furthermore one has to require a law of large numbers for the power spectrum P A which 
can be expressed in the form 

lhnW*(l n h,kulr 2 h,k2) = W(r u k 1 )W(r 2 ,k 2 ) (1.24) 

X— *0 

for r\ + r 2 . 

In general, the assumption of strict Gaussianity is too strong. Because of the anhar- 
monicity, unavoidably there will be small errors. To allow for them we call an arbitrary se- 
quence of probability measures, locally stationary, if the Gaussian property is gained 
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only in the limit A — * 0. More precisely, for our model we may assume that all odd mo- 
ments of (-) A vanish, since this property is propagated in time. For the second moments 
(a(k)a(k')) A , (a(k)*a(k')) A we form the rescaled one-particle Wigner function as in (11.201 1. 
(11.21b . The rescaled Wigner function for (a(k)a(k')) A is assumed to vanish as A — > 0, 
while the rescaled Wigner function for (a(k)*a(k')) A satisfies (11.23b . ( 11.24b is imposed 
correspondingly. In addition we require that the sequence of local measures, close to A~ x r 
[lattice units], converges to a Gaussian measure with covariance ( 11.16b . ( 11.17b and W(k) 
replaced by W(r, k). 

From ( 11.23b we infer that for a locally stationary measure it holds 

<^ha>' l = f dku){k){a(kYa{k)) A = 0(A~ 3 ) . (1.25) 

Jt 3 

The harmonic energy is extensive. On the other hand the average anharmonic potential 
satisfies ^[A{V) A - *\fA~0(A~ 3 ), which is small compared to (//ha)' 1 - 

Our real goal is to understand the time evolution in the limit of small A. We impose 
the initial measure to be Gaussian and locally stationary. As a simple, but useful first 
step we set V = and, for a short intermediate stretch only, consider the linear dynamics 
generated by H^. By linearity this dynamics preserves the Gaussian property. Since the 
initial state has a slow variation on the scale A -1 , the limit A — * for the one-particle 
Wigner function is a particular case of the semiclassical limit for a linear wave equation. 
Thus the appropriate time scale is also of order A~ l and on that scale 

hmW\[rh,k,Ar 1 t) = W(r,k,t). (1.26) 

/t-»0 

(On the left hand side A~ l t is the time in microscopic units, the units of ( 11.131 ). while on 
the right hand side t refers to kinetic times, the units for ( 11.15b ). The limit Wigner function 
is governed by the transport equation 

|-W(r, k, t) + ^V k co{k) ■ V r W(r, k,t) = 0, (1.27) 
of 2n 

which corresponds to the free motion of fictiuous particles, the phonons, with kinetic 
energy u>. Of course, ( 11.27b has to be solved with the initial condition W{r, k, 0) = W(r, k). 

We refer to |fl6l [171 for a complete coverage of the semiclassical limit for lattice 
dynamics with an arbitrary unit cell. Random initial data are studied in IfTSl , where only 
rather mild mixing conditions on the initial measure (-)' 1 are imposed. Let (■)* be the 
measure at time t as evolved according to the dynamics generated by //ha- It is proved that 
for any t > the sequence of measures (-)j_ lt is locally stationary with a Wigner function 
satisfying the transport equation ( 11.271 ). In this sense local stationarity is dynamically 
generated. Such results require, as does the kinetic limit, that phonons propagate with 
nonzero velocity. To say one has to demand that V^oi + a.s.. By our assumption, to 1 is 
real analytic on T 3 and one only has to exclude the case co = const.. Nevertheless, since 
momentum space is T 3 , there will always be submanifolds where V^w = 0, which sets an 
extra technical difficulty. 
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We return to the case of interest, namely adding the anharmonicity yfXv. Then, if 
as before the initial measure is required to be Gaussian and locally stationary, the time- 
evolved measure {-}f' A is no longer Gaussian. On the other hand, the linear dynamics does 
not tolerate deviations from Gaussianity, as demonstrated in [18|, and the strength of the 
anharmonicity is taken to be small. Thus the family (-K.il, is still locally stationary. Of 
course now W(k) depends on r and t [kinetic scale]. The precise strength of the anhar- 
monic on-site potential is adjusted such that its effect is of order 1 for the scale on which 
( 11.271 1 holds. More physically speaking, the interaction strength is such that the phonon 
mean free path is of order A [lattice units]. Thus one expects the limit dl.26l ) still to be 
valid. Only the transport equation picks up a term resulting from the collisions between 
phonons. 




Figure 1. Time evolution of a probability measure under kinetic scaling. 

We summarize the kinetic framework by means of a somewhat schematic diagram. 
The encircled set is supposed to represent the set of "all" probability measures for the 
random field a(k). We fix the interaction strength as < A «: 1. The manifold-like set 
denoted by S' 1 is the set of all locally stationary measures up to some /l-dependent scale 
of precision. The hamiltonian dynamics generates a flow on the space of measures. The 
set S' 1 is attractive, in the sense that fairly rapidly, times of order 1, the measure will be 
close to S' 1 . On S' 1 the time evolution is slow with changes on the kinetic scale 0(A~ l ) 
and the defining Wigner function is governed by a kinetic equation. For t » AT further 
dynamical phenomena set in which cannot be captured through the kinetic framework. 

2 The Boltzmann equation 

The kinetic framework provides us with a tool for an educated guess on the form of 
the collision operator. An important check will be that the resulting equation has the 
physically correct stationary solutions and satisfies an //-theorem for the total entropy, 
respectively a local semi-conservation law with a positive entropy production. 

A more demanding issue is to support the educated guess by mathematical arguments 
which establish the assumed existence of limits and local stationarity. 



8 



Herbert Spohn 



Let (-)f' A = (■); denote the measure at time t under the dynamics generated by H with 
(')o = (■) G ' A - The two-point function satisfies 

^Mpfaiq)), = i(cj(p) - oj(q))(a(pTa(q)) t + ^AF(q,p, t) . (2.1) 
at 

F is cubic in the a-field. We integrate ( 11.131 ) in time, symbolically 

a(t) = e~ i0Jt a + VI f e - iai( '- s) a(sf , (2.2) 
Jo 

and insert ( 12.21 ) in the expression for F. There are then cubic terms, proportional to y~A, 
which vanish upon averaging, and sixtic terms, proportional to A. At this point we use, 
as argued, that locally the measure is approximately Gaussian and factorize the 6-point 
function according to the Gaussian rule for moments. We also rescale space-time by A~ l 
and switch to Wigner function coordinates. In this form one can take the limit A — > 0. 
The first term on the right of ( 12.11 ) yields the free flow while yfAF results in the collision 
operator with a cubic nonlinearity. 

Details of such computations can be found in (9] \l9l . Here we only quote the result, 
which is the Boltzmann equation 

d 1 

— W(r, k, t) + —V k u{k) ■ V r W(r, k, t) = C(W(r, t))(k) . (2.3) 
ot 2n 

The collision operator is local in r, t, which our notation is supposed to indicate. C is a 
nonlinear functional of k i-> W(r, k, t) at fixed r, t. For the quartic on-site potential, C is 
cubic and defined through 

C(W)(k) = 12tt V I dkidk2dk 3 (\6ucj X (x)2a)^y l 

o-uo-2,<r 3 =±l ^ T9 

x6(k + criki + a"2^2 + crj,kj,)5{u) + ctioji + (T2OJ2 + cr-jOJi) 
x(WiW 2 W 3 + W(o-iW 2 W 3 + Wio- 2 W 3 + Wi W 2 cr 3 )) , (2.4) 

where we use the shorthand Wj = W{kj), uj = 0)(kj), j = 1, 2, 3. The term proportional 
to W is the loss term. It has no definite sign. On the other hand the gain term W1W2W3 
is always positive, ensuring that the positivity of the initial Wigner function is propagated 
in time. 

As a historical parenthesis we remark that in his seminal paper Peierls uses a very 
different reasoning. The current method goes back to the mid-fifties where Green's func- 
tion techniques were developed systematically along the lines set forth by Quantum Field 
Theory. The textbook expositions of the phonon Boltzmann equation are not so satisfac- 
tory, at least in solid state physics. For the quantized theory Fermi's golden rule is applied 
in a spatially homogeneous, discrete mode context. While this gives the correct collision 
operator, the 'why' remains obscure. Classical lattice dynamics is conceived only as a 
Ti — > limit of the quantum theory. The situation is more favorable in wave turbulence, 
where the starting point is a weakly nonlinear system of wave equations in hamiltonian 



Weakly nonlinear wave equations 



9 



form [20|. Statistical properties are accessible and the validity of Gaussian statistics is an 
important issue ||2T1|221 . 

At present, we are very far away from a complete mathematical proof of the existence 
of the limit in Eq. ( 11.26b together with the property that the limiting Wigner function is 
the solution to ( 12.31 l. ( 12.4b with initial conditions W(r, k, 0) = W(r, k). Apparently, the 
only available technique is to expand with respect to yAV into a time-dependent per- 
turbation series. Each term of the series can be represented symbolically as a Feynman 
diagram, which represents an oscillatory integral. In the limit A — > most diagrams van- 
ish (subleading diagrams) and a few do not vanish (leading diagrams). The sum over all 
leading diagrams yields the time-dependent perturbation series of the Boltzmann equa- 
tion, where the free flow is integrated and the collision term is regarded as perturbation. 
The conjectured division into leading and subleading diagrams is explained in [23 1. 

The program outlined has been carried out only fragmentarily. The first order, propor- 
tional to yA, vanishes. To second order one expects the limit 

f dse u '- s) C(e Ls W) , (2.5) 
Jo 

where (e L 'W)(r,k) = W(r- (27r)~ 1 V k a){k)t, k). While such limit could be established at the 
level of generality discussed here, even this point has not be accomplished. Closest to the 
goal comes the very careful analysis of Ho and Landau [24], who study the same problem 
for a weakly interacting Fermi liquid on Z 3 with nearest neighbor hopping and obtain the 
analogue of ( |2.5) . Benedetto et al. l|25ll26l consider a weakly interacting quantum fluid 
in R 3 with the usual kinetic energy p 2 /2m. In our context this would correspond to the 
nonlinear Schrodinger equation, m = 1, 

i—Hx, t) = ( - ±A A - + VIV*(x, t))ifr(x, t) (2.6) 
with the effective potential 

V„(*,f) = f dy§(x-y)\i(,(y,t)\ 2 . (2.7) 

§ smoothens the interaction on a microscopic scale. The role of the a-field is taken over 
by the Fourier transform if/. The Wigner function W A (t) is precisely the standard Wigner 
function for the Schrodinger equation at time A~ x t averaged over the initial Gaussian dis- 
tribution of the i/r-field. The Boltzmann equation for the nonlinear Schrodinger equation 
(12.61 > has the same overall structure as (12.31 1. (12.4b . More precisely, now k e R 3 and 

4- W(r, k,t) + k- V,W(r, k, t) = C NS (W(r, t))(k) (2.8) 
ot 

with the collision operator 

C m (W)(ki) = I2n \ dk 2 dhdhW{h - k 2 )\ 2 26(k\ +k\-k\- k 2 4 ) 

x6(ki +k 2 -h- k 4 )(W 2 W3W4 - Wi(W 2 Wi + W 2 W 4 - W^W 4 )) . (2.9) 
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The Feynman diagrams of the nonlinear Schrodinger equation are a proper subset of 
those investigated in [25]. Using their result, one concludes that in the limit A — > the 
second order expansion term is given by (12. 5t with the appropriate adjustments for L and 
C. 



3 Stationary solutions, entropy production 

For dilute gases the /-function derived from the thermal equilibrium distribution of the 
-particle system is a stationary solution of the kinetic equation. If one identifies the 
entropy with the logarithm of phase space volume associated to some given /-function on 
one-particle phase space, then the entropy functional turns out to be 

S(f) = - J rfv/(v)log/(v), (3.1) 

in units where Boltzmann's constant k% = 1. S (/) is increasing in time and constant if 
and only if / is Maxwellian (at least in the spatially homogeneous case). Thus the reader 
may wonder whether the phonon Boltzmann equation has similar properties. 

To determine the collision rule for phonons, one has to solve the conservation laws of 
energy and momentum, see Eq ( 12.4b . For pair collisions, i.e. £j=i °~j - -1 > this amounts 
to 

co(ki) + Q)(k 2 ) = o)(ki) + 0)(ki + k 2 - k 3 ) (3.2) 
and for three phonons mergers to 

a>(ki) + a>(k 2 ) + w(£ 3 ) = u)(ki +k 2 + h) , (3.3) 

where both equations implicitly define £3 = ki{k\,k 2 ). £4 is obtained from momentum 
conservation as £4 = k\ + k 2 - £3 mod 1, resp. k\ - k\ + k 2 + £3 mod 1. Mostly one has to 
work with this implicit definition and only in very exceptional cases an explicit collision 
formula is available. An instructive example is the case one space dimension. Mechanical 
particles would merely exchange labels and the Boltzmann collision operator vanishes. 
On the other hand, for phonons in one dimension with an u derived from nearest neighbor 
couplings, ( 13. 2t has a nondegenerate solution, thus providing real phonon collisions 11271 
l28l . In principle, it may happen that the energy-momentum conservation laws have no 
solution at all. Then the collision term vanishes. In case of several solutions one has to 
sum over all of them. Note that in ( 12.41 i energy conservation cannot be satisfied for the 
term with crj — 1, j — 1, 2, 3. This term had been added only to have a more symmetric 
looking expression. A further instructive example are nearest neighbor couplings, for 
which the dispersion relation reads 

3 

w(fe) = (w ( 2+2^](l-cos(27rA;' / '))) 1/2 , k = (k l ,k 2 ,k 3 ) e T 3 . (3.4) 
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Then (13.3b has no solution and collision processes where three phonons merge into one, 
and their time reversal, are forbidden. As for dilute gases, there are only number conserv- 
ing pair collisions. Such a property is stable under small perturbations of a>. It also holds 
for the nonlinear wave equation ( II . lb for which a>(k) — \k\, k eM?. 

The equilibrium measure for the harmonic part is the Gaussian Z 1 exp[-/3Hh dr ] with 
Z the normalizing partition function and £1 the inverse temperature, /3 > 0. The anharmonic 
potential can be neglected in the kinetic limit. The corresponding Wigner function, Wp, is 
given by 

W = ~R~nA ' (3 ' 5) 

in the limit of an infinitely extended lattice, compare with ( 11.16b . Using (12.41 >. it follows 
that 

C(Wf}){k) = ^f ^ I dk\dk2dkj,(oJOJ\<jj2 (jJ i)^8{oj + o~ \U)\ + o~ 2^2 + cr 3W3) 

o-i,<r 2 ,o- 3 =±l ^ T3 

x6(k + a-\k\ + cr 2 ^2 + cr^k^)^ (0)0)10)20)3) 1 (to + o-\lo\ + cr 2 W2 + o~3U>3) 
= 0, (3.6) 

as expected. The issue whether there are further stationary solutions will be discussed 
below. 

Following Boltzmann, to define the entropy one first has to identify a family of macro- 
scopic observables. If for simplicity, and in fact for the remainder of this section, we 
restrict ourselves to the spatially homogeneous situation, then in spirit the macrovariables 
are the phonon numbers a{k)*a{k), k e T 3 . To be more precise one has to employ the 
limit procedure of Boltzmann, the only difference being that here we have to take into 
account the harmonic couplings. We partition the torus T 3 into cubes A ; of side length 
6, j = 1, . . . ,M 3 , 5M = 1. We consider a finite lattice volume, { 3 , which is equivalent 
to discretizing the torus T 3 to (IV) 3 with grid spacing £~ l . The proper macrovariables are 
then 

Hj= Y a{k)*a(k), ;'=1,...,M 3 . (3.7) 

*eA,n(T,) 3 

For a given Wigner function W the cell in phase space at precision f 3 s is defined by the 
conditions 

{{q,p)e{R 1 f\{\e j -s)<H j {q,p)<e\e j + s), ; = 1, . . . , M 3 } (3.8) 

with 



1 



6 3 I dkWik). (3.9) 



The entropy S(W) of W is defined as times the logarithm of the Lebesgue measure of 
the set in ( 13.8b upon taking limits in the following order: I — > 00, s — > 0, 6 — > 0. The net 
result is 

S(W)= f dklogW(k), (3.10) 
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up to a constant independent of W. 

Before computing the rate of change of the entropy let us introduce the notion of a 
collisional invariant. This is a function <// : T 3 — > R which satisfies either one of the 
following functional equations. 

(i) (pair collisions) 

Iff(ki) + lff(k 2 ) = llf(k3)+lff(k 1 +k2-k 3 ) (3.11) 

on the set {(k\, ki, k$) s T 9 [ <x>(k\) + loiki) = <o(k$~) + cj(k\ + k2 - feO). 

(ii) (three phonons merger) 

.K&i) + i/,(k 2 ) + ifr(h) = iKh +k 2 + h) (3.12) 

on the set {(k\, ki, k$) e T 3 | a>(ki) + u(k 2 ) + u(ks) = cj(k\ + k 2 + £3)}. 

If W evolves according the spatially homogeneous kinetic equation (12.3b . then 



j ( S(W) = | (IhWa-.r'CtWHk:) 



-X' 

— Yin I dk\dk 2 dkidk^(\6(jj\CL> 2 (jjT,a)i[) l 5{ ^ o~jU)j)5( ^ 

(Ti,o"2,cr3,a"4-±l _/=l y'-l 

xWf'o-Ko-i W^^Wi + W10W3W4 + W1W20-3W4 + WiW 2 W 3 cr 4 ) 

r 4 4 

= 37r ^ I dk\dk 2 dk-}dk4(\6<jL>\(±> 2 (jjT,CL>4Y l 6( ^ o~jUj)5( ^ 

O"! ,0"2,a"l,0"4=±l ^* 7=1 7=1 

4 

xWiWaWaW^^cr/W; 1 ) 2 . (3.13) 

7=1 

We conclude that the entropy production is non-negative and that dS(W)/dt = if and 
only if 1 /W is a positive collisional invariant. 

Vice versa, if ^ is a collisional invariant and tfr(k) > 0, then setting W(k) = l/iff(k) one 

has 

C(W)(k) = 0. (3.14) 

Thus stationary solutions are uniquely characterized by 1 /W being a positive collisional 
invariant. In brackets we remark that a complete discussion would have to specify in 
which function space one searches for collisional invariants. 

Collisional invariants are studied in [29], where the following result is proved, co is 
assumed to have no fiat pieces in the sense that the set \k e T 3 1 det Hessa>(k) = 0} has at 
most the dimension two. Here Hess is the Hessian of a> as a 3 x 3 matrix. Let 1// : T 3 — > E 
be measurable, dk\ip(k)\ < 00, and a collisional invariant under pair collisions, in the 
sense that (13.111 1 holds Lebesgue almost surely. Then <// is necessarily of the form 

iff(k) = a + cu>(k) (3.15) 



with some constants a, c e R. By our previous discussion this implies that the kinetic 
equation has a two-parameter family of stationary solutions. If ( 13.151 1 is inserted in the 
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3-phonons merger (13.121) . then a — necessarily. To have thermal equilibrium as only 
stationary solutions, there must be 3-phonons mergers on a set of full dimension. A more 
physical mechanism would be to include a cubic on-site potential, which appears naturally 
in an expansion with respect to the anharmonicity, unless there are special symmetries 
which would make this order to vanish. For the cubic potential a collisional invariant 
must satisfy the 2-phonons merger condition, which is (I3.12l i upon dropping k^. Again, if 
there is a set of full dimension of 2-phonons mergers, then the only stationary solutions 
are the ones corresponding to thermal equilibrium. 

Our discussion has interesting implications on the very widely studied Fermi-Pasta- 
Ulam problem [30|, see the special issue OTl at the occasion of its 50-th anniversary. 
Fermi, Pasta, and Ulam considered chains of anharmonic oscillators and studied by a nu- 
merical iteration scheme the relaxation to equipartition, i.e. to thermal equilibrium. To 
their surprise, such an approach did not take place, at least not on the time scale of the 
computation. One explanation comes from the KAM theory which demonstrates that in 
part the invariant tori of the linear dynamics persist under small nonlinear perturbations. 
Thus a Lebesgue typical phase point may not explore the full energy shell. A second 
explanation relies on the fact that for a special choice of the nonlinearity the Hamilto- 
nian system remains integrable and admits propagating soliton solutions. Also, for the 
FPU model there are special solutions, the breathers, which may impede the relaxation to 
thermal equilibrium. 

Kinetic theory operates in a different part of phase space. The energy is proportional 
to the length of the chain and the nonlinearity is weak. In this regime the phonon picture is 
precise. Propagation is encoded by the dispersion relation co, while the relaxation comes 
from the phonon collisions. For the FPU chain one has u{k) = (1 - cos(27rfc)) I/2 , tel. 
Thus, ignoring special fc-values, phonons propagate ballistically with non-zero speed. On 
the other hand, the collision rule depends on the precise form of the nonlinearity. For 
the FPU a-chain the couplings are cubic as (qj+i - qf) 3 . In this case the collision term 
vanishes which signals poor relaxation. For the FPU yS-chain the couplings are quartic as 
(qj+i - qj) 4 . A detailed analysis shows that the energy current correlation, as based on 
kinetic theory, has a slow decay as f~ 3 ^ 5 11271 [33l . a result with good numerical support 
ll32l . If instead one chooses as nonlinearity a quartic on-site potential, as q 4 j, then the 
energy current correlation decays exponentially 11281 . In this case kinetic theory predicts a 
rapid convergence to equilibrium, as is well confirmed by molecular dynamics simulations 
ED [34 1, also away from the limiting kinetic regime. 
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